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An important problem associated with any mission involving a space- 
craft or aircraft is that of correcting path errors in such manner that 
prescribed terminal conditions which define successful completion of the mis- 
sion are satisfied. Determination of these corrections is referred to as 
guidance. Many vehicles must be controlled by automatic or semi-automatic 
devices which operate without the direct intervention of a man. For non- 
tracked vehicles, these devices must control the path and speed of the vehicle. 
It is this type of vehicle (spacecraft or aircraft) which will be considered. 
There are several definite steps in the solution of the guidance problem. 

First, a nominal trajectory must be defined. Next, the errors in a current 
state can be evaluated by comparison to the nominal trajectory. Using infor- 
mation from the nominal trajectory plus the known state errors, guidance in- 
formation can be generated to drive the terminal state toward the desired end 
conditions. 

This study compares two variations of the minimum distance guidance 
algorithm (open and closed loop forms) with a nonlinear reoptimization scheme 
which is the ideal guidance procedure. This study seeks to determine the ef- 
fects of loop closure on the performance of the guidance algorithm. For the 
numerical simulation, a minimum time low- thrust Earth-Mars transfer is used 
as an example problem. 

The author wishes to express his deep appreciation to Dr. W. T. 

Fowler for serving as supervising professor and for his helpful suggestions 
and encouragement during the course of this investigation. He also would 
like to express his appreciation to Dr. D. G. Hull for serving on the super- 
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ABSTRACT 


A comparison is made of open and closed loop applications of a 
second order guidance algorithm, using the minimum distance strategy. A non* 
linear reoptimization procedure is used as the ideal guidance history. The 
system model used for the comparison is a low-thrust vehicle performing a 
minimum time, three-dimensional, heliocentric Earth-Mars transfer. 

For the example problem considered, closed loop guidance proves to 
be much more accurate cn satisfaction of the final state than the open loop 
procedure. On the other hand, closed loop guidance proves to be much more 
vulnerable to perturbation by highly nonlinear regions in the trajectory. 
Finally, the results indicate that for this problem the best loop closure in- 
terval is at each integration step, about one day, or more often, if possible. 
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[JET or EYMBOLE 


Tiir' following list tubulates all significant symbols used in the 
main text. Each symbol is accompanied by a brief description and the equa- 
tion number where the symbol is introduced. 


MATRICES : 

The size of a matrix is ii. Heated in the statement immediately 
following the symbol. The following specific indices are used. 

n - number of state variables 
m - number of control variables 

n - number of initially specified constraint relations 
1 - number of terminal constraint relations 


[A] n*n matrix of partial derivatives 

(X(t,t^)] 2n*(n+q) matrix of partial derivatives resulting from the n+q 

backward integrations of the 2n vector of perturbation equations 

[X 1 (t,t f )J n*n submatrix of [X(t,t^)J 

[X 2 (t,t^)J nxq submatrix of [X(t,t f )] 

t x 3 <t»t f )J n*n sub matrix of [x(t,t^)J 

lX^(t,t^)J nxq submatrix of tX(t,t^)] 

■xn matrix of gain functions 
mxq matrix of gain functions 
mxn matrix of gain functions 
mxn matrix of gain functions 


ca 2 i 

[a 3 j 


[\1 

[A(t|t)I 2n*2n matrix of partial darivatiwss resulting from tbs 2n forward 
integrations of tbs 2n rooter of parturt>ation 




nu 
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VECTORS : 

All vectors are column vectors unless otherwise noted. The vector 
size is indicated in the statement immediately following the symbol. 

f n-vector of state variable derivatives 

F 2n-vector of state variable and Euler-Lagrange variable derivatives 
h (ntl)-vector of terminal constraint relations 
M q-vector of specified terminal constraint relations 
u m- vector of control variables 
x n-vector of state variables 

z 2n-vector of state variables and Euler-Lagrange variables 
X n-vector of time dependent Lagrange multipliers 
v q-vector of constant Lagrange multipliers 
p (n-q-i l)-vector of terminal constraints 

SCALARS : 

c propellent exhaust speed relative to vehicle 
H variational Hamiltonian 
I auxiliary functional to be extreaised 
J performance index 

K terminal constraint obtained from variational analysis 
■ instantaneous mass of vehicle 
■ e initial mass of vehicle 
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I' .juxilinry I utict ionfll oqu-il to ij> + v M 

r imt.jntaneou:; radial distance of vehicle from origin of heliocentric 

coordinate system 

t independent variable, time 

T thrust of the ion engine 

6 constant mass ejection rate 

u gravitational constant for the Sun 

<f> performance index 


SUPERSCRIPTS : 

< ) differentiation with respect to time 

T 

( ) matrix transpose 

( ) -1 matrix inverse 

( )* indicates evaluation on the reference trajectory 



SUBSCRIPTS : 

( ) evaluation at the initial time 
o 

( )f evaluation at the terminal time 



( evaluation at some time while the mission is in progress 
( )„ Mars 

( ) E Earth 

ABBREVIATIONS : 

AU astronomical unit 

HOG minimum- distance guidance 
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1. INTRODUCTION 


1.1. Preliminary Description 

As man extends his horizons to include exploration of the solar sys- 
tem, large mission durations are encountered which would present some prob- 
lems for a manned spacecraft. A logical step would be to shorten these mis- 
sion durations (transfer times) with the vise of some additional thrust during 
the transfer. The most likely candidate for this additional thrust is a low- 
thrust, continuous operation engine. This might be an ion or plasma-type 
engine. The force produced by this thrust requires the development of a 
method to control the direction of thrust so that the terminal conditions of 
the mission will be satisfied. This type of mission would be characterized by 
four phases; a) ascent and escape phase, b) midcourse phase, c) approach 
phase and d) terminal maneuver phase. 

The midcourse phase is a heliocentric transfer with the spacecraft 
engine thrusting continuously at some low level. For the class of trajector- 
ies that does not contain coast periods, it is desirable to minimize the 
trans er time. This will allow a minimum amount of fuel to be carried through 
the ascent and escape phase. The constant magnitude thrust can be controlled 
in direction so that the transfer time is a minimum. The purpose of this 
study is not to look at different trajectories but to compare the performance 
of two guidance algorithms on the same trajectory. A nominal minimum time 
transfer trajectory and the associated control can be generated and stored 
before the comparison is started. This nominal trajectory was generated by 
Hart (1) and is covered in detail there. 

1.2* Problem Statement 


Th# fu!4fcO» 
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be thought of as 





thftt of dtlifilnlTH thi 




set of control changes necessary to get the spacecraft be" '•o the nominal 
path after an off-course or off-nominal condition is r» - ognized. Hart (1) 
and others indicate that rather than trying to get back on the nomrnal path, 
one should find a new optimum path to satisfy the mission arm conditions. 

In the absence of an onboard computer capable of carrying out a total trajec- 
tory reoptimization, one might use a guidance routine. There are several al- 
gorithms under consideration for this type of problem. All of these methods 
are derived from the Calculus of Variations. A general comparison of several 
algorithms which might be used for a reoptimization approach was made by 
LewalJen (5). The guidance algorithm considered in this study was presented 
by Hart (1), along with several others for comparison. In Hart's study, the 
minimus distance guidance algorithm was found to be the best among those he 
considered. To the author's knowledge, all of the previous work done with 
the minimus distance guidance algorithm has been carried out employing the 
algorithm in an open loop form. This study applies the minimum distance 
guidance technique in both open and closed loop forms and compares these re- 
sults. The basis for comparison is the result of the nonlinear reoptimiza- 
tion procedure used by Hart (1). The state history generated by the minimal 
distance algorithm and that generated by the nonlinear reoptimisation pro- 
cedure are compared to try to determine the performance of the cloeed loop 
versus the open loop application. 

Several problems enc oun tered in the application of the method ere 
dis cu ss e d. The comparison of the open sod cloe ed loop analyses consists of 
calculations based on perturbations in position and velocity components at 
samuel points along the tr a jectory. At each of these points, di ffer e nt loop 
closure period ate co mp ar ed, from this information, seme optimum value for 
the loop closure period can be determined. As a ssoondary goal, if the elrmed 
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loop analysis proves to bo hotter, then some optimum value for the loop clo- 
sure period will be sought. 


1.3. Notation Convention 

All vectors are assumed to be column vectors and will be denoted 
by a bar; i.e., x , z . All matrices will be denoted by brackets; i.e., 
[A] , [H--] * The su P erscr iP t T indicates the transpose of a vector or 

matrix, and the superscript -1 indicates the inverse of a square matrix. 
All first partial derivatives of scalars with respect to a vector are now 
vectors and are denoted as : 


9H u / 3H 3H 9H \ 

7 * H ; = 57iT •" 3x~ I 

2x \ 1 2 n I 

All second partial derivatives of scalars with respec* ctors are mat- 

rices , i.e. , 



All first partial derivatives of vectors with respect to vectors are matrices 


formed in the following manner: 




if 

1 

... _I 


*n 

• 

• 

• 

• 

if 

3f 


... — &. 

*1 

*n 


The variation of x is denoted by 6x and the total changa in x is denoted 
by 6 k • In general. Ax * 6x * mdt where x ii 



?. MATHEMATICAL MODEL 


The nominal trajectory employed in this study is a minimum time, 
low-thrust heliocentric transfer from Earth to Mars. The equations of motion 
describing the trajec *ry are expressed in terms of a heliocentric rectangu- 
lar coordinate system, x 4 , x^ , x^ . The x^ axis coincides with the 
line of the ascending node of Mars, the x^ axis lies in the ecliptic plane 
and the x^ axis coincides with the angular momentum vector of the Earth 
with respect to the Sun. For simplicity, the spacecraft is given the velo- 
city of the Earth at the beginning of ; mission rather than escape velocity 
relative to the Earth. Gravitational effects of the Earth are then neglected. 
Details of the arrival at Mars will be determined by specific mission require- 
ments. For the problem under consideration, the end conditions are taken to 
be the position and velocity of a point on the orbit of Mars at the final 
time. 

Gravitational effects of Mars, along with perturbing effects of the 
other planets, have been neglected although these forces would have to be in- 
cluded in an actual mission analysis. These forces are neglected so that the 
effects of the guidance method can be studied without complications due to 
the effects of small perturbing forces in the system. 

The orbit of Mars is assumed to be an ellipse with an eccentricity 
of e « .093393 , a semi-major axis of a * 1.523691 AU (Astronomical 
Units) , and an angle of Inclination of 1 * *0??289 radians . The argu- 

ment of perihelion la taken to be a. 992927 radians , and for the coordinate 
system empl oyed, the argument of the seconding node is 0.0 . This model is 
the sane as that used by Fovler (3), Hart (1) and others. 
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3. GUIDANCE ALGORITHM 


3.1. Second Order Guidance Procedures 

The minimum distance guidance algorithm and the perturbation func- 
tions guidance algorithm are both based on the second variation of the Cal- 
culus of Variations. Both procedures calculate the same guidance information 
The difference between the two procedures is the point at which the infor- 
mation is applied on the trajectory. Hart (1) gives a more detailed treat- 
ment of several second order algorithms. 

The differential equations describing the system to be controlled 

are 

x = ?(x,u,t) (3.1.1) 


where 

x is an n-vector of state variables, 
u is an m- vector of control variables, 
t is the scalar independent variable, time, and 
f is an n-vector of known functions. 

Assume that a nominal trajectory has been obtained which minimizes a perfor- 
mance index 


J = *(x(t f ),t f ) (3.1.2) 

and satisfies the terminal constraints 

M(x(t f ),t f ) * 0 (3.1.3) 

where 4 is a scalar quantity evaluated at the terminal time and M is a 
q- vector of terminal constraints. 

This trajectory oust satisfy certain differential equations; these are 
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X = f(x,u,t) 

(3.1.4) 

’ T 

A = -H- 

(3.1.5) 

X 

H- = 0 

(3.1.6) 

u 


-Tri 

is an n-vector of multipliers, and H 

= A f 


For t < t < tj. where 
Formulation and solution of this problem is based on a Calculus of Variations 
approach. Since x(t Q ) and t Q are known. Equations (3.1.4) through (3.1.6) 
and subject to the following boundary conditions: 

At t 


x(t ) = x 

o o 


(3.1.7) 


At t. 


M(x(t f ),t f ) 


x- T (t f ) 


(A T f + * t + v T M t ) t 


= 0 


+ "Vt. 


= 0 


(3.1.8) 

(3.1.9) 

(3.1.10) 


After considering a perturbation of the nominal trajectory due to changes in 
initial and/or final conditions, Equations (3.1.4) through (3.1.6) become 

5t (6i) = + [H \u ]6 * * (3.1.11) 

(«X) * -[H~J6x - [H--]6u - [H-rl«X (3.1.12) 

[H~J«S ♦ ♦ [HjjHX ■ 0 . (3.1.13) 

The boundary conditions. Equations (3.1.7) through (3.1.10), become 


dx(t-) » dx 
o o 


(3.1.14) 


:r : 
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dM f = {[M-]6x + [M-]fdt + M^dt}^ = 0 


(3.1.15) 


•' .. 4 
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A total variation of M(x(t^),t^.) rather than a partial variation must be 
considered since the final time, t^ , as well as x(t^) is allowed to change. 
Thus , 

6x^. = { [P — ]6x + ( [P — ]f + (P -) T + H-)dt + [M-] T dv} (3.1.16) 
f xx J xx tx x x J t f 

{ 6 A T f + (H- + P -)6x + (H + P -f + P )dt + dv T M } = 0 (3.1.17) 

x tx t tx tt t t ^£- 

-T- 

where P = <j> + v M , 

This set of Equations, (3.1.11) through (3.1.17), must be satisfied over the 
entire trajectory, for t^< t, < t^ . 

Now, assuming that [H--] is nonsingular. Equation (3.1.13) can be 


solved for 6u to give 


Su ^ -[H— ] -1 { [H — ]6x + [H- t ]6A> 

1 uu J 1 ux J 1 uA J 


(3.1.18) 


Using this result, Su can be eliminated from Equations (3.1.11) and (3.1.12), 
This will yield a set of 2n first order linear differential equations, with 
time dependent coefficients in 6x and 6 A . For simplicity, let this sys- 
tem of 2n equations be written as 


4- * 

dt U 


[$!•;•$!] CM 


(3.1.19) 


where 


[A,(t)l - [Hjjl - IHjGlWsGl'hHQjl . 

»,<«>) • -(HjjHBjjr'lHjjl . 
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[A 3 (< )l 

- 't'fei 

+ [H — ] [II- - ] 
1 XU uu 

‘l“ux’ 

[A 4 (t)] 

= -‘"H 1 

+ [H — ] [H — ]“ 
1 XU uu 
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Numerical integration of Equation (3.1.19) with proper boundary con- 
ditions will now determine 5x and 6 A . These results can be used in Equa- 
tion (3.1.18) to evaluate the control variation history, 6u(t) . The bound- 
ary conditions can be obtained from Equations (3.1.15) through (3.1.17). The 
n+q+1 equations contain the 2n+q+l unknowns 6x^ , <5A^ , Av and At^ . 

If the n+q+1 equations are linearly independent, then there are n indepen- 
dent quantities among the unknowns. Therefore, if n of the unknowns are 
chosen arbitrarily and independent of one another, the 2n+q+l unknowns are 
unique determined. Choose the q components of Av and n-q of the compo- 
nents of dx^ as the arbitrarily chosen unknowns. By specifying each of 
these values, all 2n+l components of 6x^ , 6X^ and At^. can be determined. 
This now allows backward integration of Equation (3.1.19) from t^ to t^ . 
The value of t^ may be any arbitrary value in the interval t Q < t^ < t^ . 

The n arbitrarily chosen unknowns and the q components of AM^. are 
all specified quantities at t^ . Combining these vectors will give an n+q 
vector of specified quantities* Equation (3.1.19) is now integrated n+q 
times with n+q linearly independent starting conditions. The n+q vectors 
of solutions for 6x(t) and 5A(t) form a (2n*n+q) transition matrix. 

After some manipulation the desired solution for 6x(t) and 6A(t) is ob- 
tained. 
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where 


[x^(t,t^-)] is an n-n matrix, 


[x (t,t^-)] is an n x q matrix, 

[x^(t,t_p)] is an n*n matrix, 

[x^(t,t^)] is an n*q matrix, 

0 is a q*n null matrix, and 
q*n 

1 is a qxq identity matrix. 

qx q 

When Equation (3.1.20) is substituted into Equation (3.1.18), the control var- 
iation, 6u(t) , in terms of the state error, 6x(t^) , and terminal constraint 
error, AM^ is obtained. Making this substitution and rewriting gives the 
following guidance equation 


6u(t) 


[A 1 (t),A 2 (t)] 


[x 1 (t 1 ,t f )] _1 


( 3 . 1 . 21 ) 


[x 1 (ti,t f )]“ 1 [x 2 (t 1 ,t f )] 


I’.i 


where 


[q(t ) ] = [HQQ]" 1 {[H”][x 1 (t,t f )] + [H-^][x 3 (t,t f )]} 


[A 2 (t)J = [H- 5 ]-l{[H~][x 2 (t.t f >] + [H a j][x 4 (t,t f )3> 


After the nominal path has been determined which minimizes 
J s ♦(x(t-),tj) and satisfies the necessary boundary conditions, the com- 
putational procedure for this algorithm can be sunnier! zed as follows: 
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Solve the n+q+i Equations (3.1.15) through (3.1.17) n+q 
times. From the solutions of these equations, n+q initial 
conditions are obtained for the integration of Equation 
(3.1.19). 


2. Now integrate Equation (3.1.19) n+q times from t^. to 
tj . Store the integrated values at each step. 

3. The gain matrices, [A^(t)] and [ ( t ) ] can now be eval- 
uated. Substitute these results into Equation (3.1.21) to 
give the control variation history, 6u(l) . 

This presents a general algorithm for generating second order guid- 
ance information. The perturbation functions approach and the minimum dis- 
tance procedure differ only in the way this guidance information is applied. 


3.2. Variations in the Procedures 

The state variations for the perturbation functions guidance al- 
gorithm and the minimum distance algorithm are determined in the same manner. 
For the control variation, however, the calculation is the same, but the guid- 
ance information is applied in a slightly different way. 

In the minimum distance guidance algorithm, an index time is used 
for calculating the gain matrices which minimizes a metric function of the 
difference between the current perturbed states and the states on the nomi- 
nal trajectory. In other words, the index time is taken to be the time of 
the closest point on the nominal trajectory. The perturbation functions 
guidance algorithm uses the time based on the time of the perturbation on 
the nominal path for calculation of the gain matrices. This procedure was 
originally proposed by Powers (2) and is covered in detail in Hart (1). The 
correction to the time for calculation of the gain matrices is given as 


v**ffV'-**’ , i 
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’’ * .♦. \ 

-x'. (x. - x? ) 

6 ill 

At = ? <3—^ — 

L = i( x. (x. - x. ) - x. - ' 
aii 1 

c c cj 


( 3 . 2 . 22 ) 


where 


( ) indicates evaluation at the current time and 
c 

( indicates evaluation on the nominal path. 

This calculation is based on the position only, as indicated by the summation 
from four to six. A value of At is calculated at each point along the nom- 
inal path beyond t^ . The index time, whici. is used for calculating the 
gain matrices, is now determined by t + At . The nominal state variables 
at this index time are used to calculate the gain matrices. In order to make 
an evaluation of the performance of the algorithm, the guidance information 
obtained from the technique is added to the nominal control program. The 
nonlinear equations of motion are then integrated, using this augmented con- 
trol history as the control program. This integration gives the state his- 


3.3. Closed Loop Application of the Algorithm 

The previous section covers only open loop control, using the min- 
imum distance guidance procedure. A closed loop modification of this algo- 
rithm is now considered. The procedure remains basically the same as far 
as solution of the n+q+1 Equations (3.1.15) through (3.1.17) and the back- 
ward integration of Equation (3.1.19) from t f to t^ . The gain matrices, 
(A 1 (t)J and [A^(t) J , are now evaluated. This allows calculation of the 
control variation history and the state variation history. The guidance 
deviatim* predicted at t^ are now added to the nominal guidance history. 
The nonlinear equations of motion for the vehicle are integrated using this 
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.-lugmrntod guidance i nf ormat i on . At thir. point, rather thin allowing the guid- 
ance information to bo used throughout the root of t ho mi ;;r. i on , one can cre- 
ate a closed loop system to update the guidance. The guidance loop is closed 
at periodic intervals by evaluating the state errors at this time. The gain 
matrices are then re-evaluated. A new control variation history is calcu- 
lated, along with a new state deviation history. The new augmented control 
history is formed, and the nonlinear equations of motion are once again in- 
tegrated. The advantage of this closed loop method is that it continually 
updates the control variation history as the mission proceeds. The computa- 
tional procedure for the minimum distance algorithm, open and closed loop 
can be summarized as follows : 

1. Solve the n+q+1 Equations (3.1.15) through (3.1.17), n+q 
times. From the solutions of these equations, n+q initial 
conditions are obtained for the integration of Equations 

(3.1.19) . 

2. Now integrate Equation (3.1.19) n+q times from t^ to t^ . 

Store the integrated values at each step. 

3. Evaluate the gain matrices, ( ( t ) ] and [ A^(t ) ] , using 
the states for the closest point on the nominal trajectory. 

4. Substitute (Aj(t)] and [A 2 (t)] into Equation (3.1.21) to 
give the control variation history, 6u(t) . 

5. Evaluate the state variation history, Ax(t) from Equation 

(3.1.20) . 

6. Form the augmented control history and integrate the nonlinear 
equations of motion using the new control. 

7. a. If the control is to be open loop, the algorithm is fin- 
ished. 

b. If the control is to be closed loop, integrate to the 
tins of the next loop closure, evaluate the state deviations 
and go back to Step 3. 


NOMINAL PATH 

PERTURBED PATH: 

OPEN LOOP 

CLOSED LOOP 
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FIGURE 3. 

TIME PARAMETERS ASSOCIATED WITH 

MINIMUM DISTANCE GUIDANCE, 
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algorithm were compared. The performance of each method was referenced to 
the ideal guidance history produced by the nonlinear guidance scheme. The 
results of this procedure are considered to be the new nominal or ideal guid- 
ance history considering the known state errors. By integration of the non- 
linear equations of motion, position and velocity histories are obtained. 

The Euclidean norms of the differences in the position and velocity compo- 
nents were computed at each point along the three trajectories generated by 
these methods. The deviations in velocity and position between the minimum 
distance techniques and the ideal guidance were then plotted as a function 
of time. The open loop form of the minimum distance technique was found to 
perform as indicated by Hart (1). However, it is of interest to discuss some 
of the problems encountered with the implementation of the closed loop control 
algorithm. One problem area was the period at about 100 days into the mis- 
sion when it is necessary to turn the spacecraft around for the deceleration 
for Mars rendezvous. The perturbed control causes the time of the turn-around 
to shift slightly. This time shift, due to the rapid turn-around causes pre- 
dictions of large control changes, as can be seen in Figure 4. These large 
control changes exceed linearity and cause the algorithm to break down. To 
avoid this problem, nominal control is used during the turn-around period. 

The second problem area occurs when the spacecraft is nearing the 
final time. At this point in the trajectory, the effect of a control change 
is small. The algorithm then predicts a large control change to cause the 
necessary results. This region is sensed by the increasing magnitude of the 
control changes. When these deviations reach a magnitude of .15 radians. 
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nominal control was used since caJculat ions of control could rut be made 
within the linear range. This limit of .IS radians was found by trying sev- 
eral values. Beyond the time at which a .IS radian control change was <ensed, 
predicted control changes quickly approached a value of abou w 3 radians. 

This limit is the same as that used by Hart (1). 

3.5. Numerical Methods 

Accuracy, time and storage requirements for these methods are all 
highly dependent on the numerical processes invo^. 'ed. Matrix inversion and 
numerical integration account for the majority of round-off and truncation 
errors. The numerical integration routine used in this study is a fourth- 
order Adams -Bahs forth predictor and a -> Ifth-order Adams-Moulcon corrector. 

A fourth-order Runge-Kutta routine is used as a starter to ontain the initial 
and three succeeding derivatives. The dependent variables are carried in 
double precision during integration in an attempt to control round-off errors. 
Matrix inversion is performed by Gaussian elimination. All numerical calcu- 
lations were carried out on the CDC 6600/6400 digital computer in the i*0RTRAN 
IV language. Storage requirements for the minimum distance algorithm is 
167,000g words while the reoptimization (ideal guidance) requires 60,000g 
words. Computation time for the minimum distance guidance was around 12 sec- 
onds. The reoptimization (ideal guidance) requires about 15 seconds *? com- 
putation time. 
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Thf iivrA cane considered involved a state error of 5 y 10“ f) AU/Day 
in the velocity component. Position deviation as a function of time is 

shown in Figure 5, for open loop and closed loop control. The loop closure 
periods are 30, 10, 5 and 1 day. A loop closure period of 30 days results 
in a larger error at the final time than does the open loop guidance. Steady 
improvement is shown, until at a loop closure period of 1 day, the minimum 
error in the final position for the closed loop method is generated. For 
simplicity, none of the plots in this study show complete results below a 
time of 80 days. The reason for doing this is that the loop closure period 
has very little effect on velocity or position deviations from initiation 
until this time is reached. A representative curve is shown in each case for 
the general closed loop results, since all of these results look basically 
the same. 

Figure 6 gives a summary of velocity deviations between the reop- 
timized trajectory and the closed and open loop techniques. The initial state 
error is again in the (velocity) component and loop closure periods of 

30, 10, 5 and 1 day are considered. Once again, the 30-day loop closure 
period is actually worse than the open loop control with gradual improvement 
as the loop closure interval decreases. A loop closure period of 1 day gives 
the least error in the final velocity deviation. 

Position deviations between open and closed loop control and the 
reoptimized for a position error of 5*1(T 6 AU are shown in Figure 7. This 
error occurs in the component. The loop closure periods considered are 
10* 5 and 1 days. Results in this case are the same as for previous cases in 
that the trend is toward increased accuracy with the smaller loop closure 
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J’or nil U\ r ;t cases i nvost igu I ed , t ho rvr.ultn are comparable. Thr 
closed loop results show continued improvement over the open loop down to the 
minimum loop closure period of 1 day. This limit of 1 day corresponds to 
the length of one integration step. 

Figures 9 and 10 show results for an en»or introduced at 65 days 
into the mission. The closed loop results appear somewhat better than open 
loop, with the errors in positio\ corrected more accurately than those ir 
velocity. 

Figure 11 gives a listing of the actual terminal errors for the 
analyses shown in Figures 5 through 8. There is actually only a small dif- 
ference between open and closed loop results at this point. Figure 12 com- 
pares the computation time for open and closed loop analyses. There is very 
little difference between the two procedures. This, in effect, means that 
the closed loop procedure is comparable to the open loop procedure as far as 
computer costs are concerned. 
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Differences Between Final States cf 
Reoptimized Trajectory and Minimum Distance Technique 


Perturbation 

Error 

Position (AU) Velocity (AU/Day) 

Loop Closure 
Period 
( Days ) 

fiXjCO) = 5.0*10" 6 AU 

1.68543 86 101x10 “ 14 

1.3356785238 x 10” 3 

Open Loop 

M 

7. 5959967022 x 10“ 5 

7. 55 81 577845 x 10“ 3 

1 

ft 

7. 59 5996 7022x10” 5 

7. 5581577845x10" 3 

5 

ft 

1.3712271024xl0” 4 

1.1165335824X10” 4 

10 

6^(0) = 5.0*1(T 6 AU/Day 

1. 6388290871x10” 4 

1. 262605942 7x10” 2 

Open Loop 


7 . 136 81 5 515 8x10” 5 

6.8088261401xl0 -3 

1 


1.0093910324xio” 4 

8. 7410210064x10" 3 

5 


1.3037531342X10” 4 

1.0469032389X10” 2 

10 


1.6468665832X10* 4 

1.2406315644X10 -2 

30 


Figure 11. ERRORS AT FINAL TIME FOR 
ERROR INTRODUCED AT INITIATION OF THE MISSION 



Velocity Perturbation in 


Component at 


Initial Point of Trajectory: 


Loop Closure 
Period 
(Days) 

T Ime 

(Sec. ) 

Open loop 

11.763 

i 

12.652 

5 

11.980 

10 

11.308 

30 

11.490 


6x 1 (0) = 5.0xl0" 6 AU/Day. 


Position Perturbation in Component at 

Initial Point of Trajectory: 


Loop Closure 
Period 
(Days) 

Time 

(Sec.) 

Open loop 

11.545 

1 

12.171 

5 

12.143 

10 

12.041 


tfx 4 (0) ■ 5.0x10”* AO. 


Figure 12. CONPOTAtlOH TINS REQUIRED 
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5. CONCLUSIONS AND RECOMMENDATIONS 


In this study, a comparison of open and closed loop applications of 
the minimum distance guidance algorithm is presented. The comparison is based 
on the differences between the results of these algorithms and the results 
generated by a nonlinear reoptimization algorithm. The trajectory used is 
for a low-thrust spacecraft on a three-dimensional, heliocentric Earth-Mars 
transfer. 

5.1. Summary of Conclusions 

The major conclusions resulting from this study are: 

1) The errors in the final state are somewhat smaller for the 
closed loop procedure than for the open loop. However, the difference is not 
large enough to indicate a choice of the closed loop procedure over the open 
loop procedure. 

2) The position errors are somewhat smaller for both procedures 
than are the velocity errors at the final time. 

3) The errors in the final state seem to become smaller as the 
loop closure interval becomes smaller. The results indicate that the choice 
of a continuous loop closure would decrease the errors in the final state but 
would increase computation time. 

4) Reoptimization appears to be the best method of guidance. Con- 
sidering the fact that the terminal errors are still relatively large* the 
closed loop procedure has only small merit above the open loop analysis. 

Since the computation times are on the same order for the 'Minimum distance 
guidance and the reoptimization (ideal guidance)* little is gained by using 
the second order guldanoe method. This result is the same as that obtained 
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by li.irl (1) in hi:; study. 

5.2. Recommendations for Further Study 

It is suggested that the study of low- thrust guidance techniques be 
extended as follows : 

1) The results of this study indicate that the nonlinear reoptimi- 
zation technique is much more accurate with only a small difference in compu- 
tation costs. Therefore, it is recommended that further effort be directed 
at perfecting the reoptimization techniques for use in low-thrust guidance. 

2) An extension of this study should be made to trajectories which 
contain coast arcs. 

In addition, trajectories for which reoptimization is not feasible 
are of interest. An example of this would be the reentry problem. 
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APPENDIX A 


Earth Orbitci. Data 

Semi-major axis, 

Eccentricity, 

Argument of perihelion, co^ 

Angle of inclination, i^ 

Argument of ascending node , 

Time of perihelion 
Period 


1.0 AU 
.016726 
0 . 0 ° 

0 . 0 ° 

n E o.o° 

Jan. 3.022307069, 1950 
365.198084 days 


Mars Orbital Data 

Semi-major axis, a 

m 

Eccentricity, « m 

Argument of perihelion, 

Angle of inclination, 

Argument of ascending node. 

Time of perihelion 
Period 


1.523691 AU 

.093393 

286.07366° 

1.84991° 

Q 0.0° 

HI 

March 17.490627, 1959 
686.868866 days 


Vehicle Thrust Constants 

Hass flow rets, 0 .00108 initial vahicla mass /day 

Exhaust speed, c .045865 AU/day 

Solar gravitational con sta nt, u 2. 96007536*10" 4 AUVday 2 
Initial ap aoecraft mass 1.0 
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Tim'* f 'i ml.,, Jon iniii'ifton K':00 noon, J.ip. '* , 

Mission duration 196 . 76694769 days 

The spacecraft leaves the vicinity of the Earth with po5;ition and 
velocity equal to that of Earth. The mission is carried out under the in- 

* fluence of the inverse gravi tational force of the Sun. After a 19^-dey helio- 

f 

! centric transfer, the spacecraft matches position and velocity with Mars. 
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APPENDIX B 

THE DIFFERENTIAL EQUATIONS OF MOT T ON AND 
EULER-LAGRANGE DIFFERENTIAL EQUATIONS 

UX 4 + T 

X. = - —T— + — COS u. cos u 0 
1 r i M 1 2 


MX 5 T 

-~x~ + TT COS u. sxn 
r 3 M 1 2 


yX 6 T . 

+ M Sin U 1 


X 4 = X 1 


X 5 = X 2 


X 6 X 3 


X 1 = 


A. = -A. 


A„ = -A, 


WA 3yx 

r 3 p“ (A 1 X 4 + A 2 X 5 + A 3 X 6 ) 


uA 3 px 

P P- U 1 X 4 + Vb + A 3 X 6 ) 


X 6 * 


yX 3 3yx 6 


r» (X 1 X 4 + X 2 X 5 + A 3 X 6 ) 


The augmented vector, z * 


(U • u 


initially 
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z - -1 .4835073*10"^ 

z 2 = 9. 2714508x10" 3 

Z 3 = °*° 

z t = 5 . 199345 x 10" 1 

z 5 = 8.3463802X10" 1 

z c = 0.0 

o 

z ? = 1.0058717029X10 1 

z g = -2.1350450338X10 1 

z g = -6. 7014133502 x 10" 1 

z 1Q = -5. 1681265185 x 10 " 2 
Z 11 = ”**• 3276 807729 x 10" 1 


-1. 3232925942 x 10~ 3 
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